Trapping polarization of light in nonlinear optical fibers: An ideal Raman polarizer 
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The main subject of this contribution is the all-optical control over the state of polarization (SOP) 
of light, understood as the control over the SOP of a signal beam by the SOP of a pump beam. We 
will show how the possibility of such control arises naturally from a vectorial study of pump-probe 
Raman interactions in optical fibers. Most studies on the Raman effect in optical fibers assume a 
scalar model, which is only valid for high-PMD fibers (here, PMD stands for the polarization-mode 
dispersion). Modern technology enables manufacturing of low-PMD fibers, the description of which 
requires a full vectorial model. Within this model we gain full control over the SOP of the signal 
beam. In particular we show how the signal SOP is pulled towards and trapped by the pump 
SOP. The isotropic symmetry of the fiber is broken by the presence of the polarized pump. This 
trapping effect is used in experiments for the design of new nonlinear optical devices named Raman 
polarizers. Along with the property of improved signal amplification, these devices transform an 
arbitrary input SOP of the signal beam into one and the same SOP towards the output end. This 
output SOP is fully controlled by the SOP of the pump beam. We overview the sate-of-the-art of 
the subject and introduce the notion of an "ideal Raman polarizer" . 



I. INTRODUCTION 

Over the past few years, the possibility of utilizing 
the Raman effect in optical waveguides as the basis for 
the development of non-linear polarizers has opened the 
way to an interesting range of potential applications, 
such as multi-channel repolarization in optical fibers, en- 
hanced amplification and even the possibility of develop- 
ing silicon-based Raman polarizers [lHEl- 

Raman-based polarization attraction falls into a broad 
class of potentially game-changing effects related to light- 
by-light control in optical waveguides. Models for such 
control are essentially nonlinear and usually imply the 
use of a high-intensity beam to modify the properties 
of the medium (for instance its refractive index or ab- 
sorption coefficient) such that propagation of a weaker 
probe beam through the nonlinearly modified medium is 
affected in a substantial and controllable way. The possi- 
bility of achieving nonlinear polarization control is rooted 
in soliton theory, namely in conservative structures such 
as the polarization domain wall solitons |6l-lll| . However 
conclusions extracted from soliton theories involving a 
medium of infinite extension can be misleading for coun- 
terpropagating waves in a medium of finite length. In 
this case the presence of boundary conditions may lead 
to solitons with a finite lifetime [12|- In such situation, 
other so-called polarization attractors representing the 
unique distribution of SOPs of the two beams inside the 
medium play a key role in the process of trapping polar- 
ization of light US- 
Different mechanisms such as photorefractive two- 
beam coupling [lij or Kerr nonlinearity [ToL Il5j have, 



over the years, proven to be capable of producing non- 
linear polarization attraction. In their initial demonstra- 
tions, all of these methods were subject to limitations 
in their application in telecommunication links: their 
response time, in the case of photorefractive materi- 
als, or the requirement of extremely high beam powers. 
Only recently results of practical relevance have emerged, 
with non-conservative schemes based on stimulated Ra- 
man [16| or Brillouin scattering [171 ]. as well as the first 
low-power lossless polarizer, consisting of a 20 km ran- 
domly weakly birefringent fiber pumped by an incoherent 
counter-propagating beam [l8[ . 

As mentioned above, here we will focus on the particu- 
lar and very promising case of Raman polarizers, in which 
the pump and signal beam propagate through a Raman- 
active medium. By way of interacting with this medium, 
the pump beam induces a phonon-mediated gain for a 
frequency down-shifted (Stokes) signal beam. The signal 
beam, co- or counter-propagating with the pump beam, 
is then gradually amplified. This amplification mecha- 
nism lies at the heart of Raman amplifiers. One degree of 
control exerted by the pump beam over the signal beam 
is the total gain experienced by the signal from input to 
the output. This degree of control is well studied in lit- 
erature and widely used in practice. Much less known 
is another degree of control - over the state of polar- 
ization (SOP) of the signal beam. The main subject of 
this study are polarization- sensitive Raman amplifiers, 
in which polarization-dependent gain (PDG), an intrin- 
sic characteristic of the Raman effect which is usually 
considered an undesirable feature in amplification appli- 
cations, can be turned into an advantage by selectively 
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amplifying only one polarization mode of the input beam. 

Signal and pump fields considered in this study are 
continuous waves (CW) or relatively long pulses, such 
that the response of the Raman-active medium is vir- 
tually instantaneous, and as such it is described by 
the instantaneous dissipative cubic nonlinearity. Mostly, 
our theoretical study is developed for silica single-mode 
fibers, though extensions to other Raman-active media, 
such as silicon are also possible Q . We shall demonstrate 
how polarization-sensitive Raman amplifiers operate in 
the regime of Raman polarizers. These Raman polarizers 
are devices that along with the function of amplification 
of light, also re-polarize the beam: the SOP of the out- 
coming signal beam is defined by the SOP of the pump 
beam, independently of what SOP the signal beam had 
at the input. In other words, the signal SOP is attracted 
(trapped) by the pump SOP. By changing the polariza- 
tion of the pump we thereby change the signal SOP. In 
this way we exercise an all-optical control over the signal 
SOP. 

In this chapter we will present the theory of Raman po- 
larizers with an emphasis on randomly birefringent fibers, 
such as the ones used in the telecom industry. We shall 
identify the conditions that are necessary for a traditional 
Raman amplifier to function as Raman polarizer, and 
characterize its performance. 



II. MODEL 

In short, we shall consider the simultaneous propaga- 
tion of two beams in a Raman-active medium. In our case 
the Raman active medium is a few kilometers long span 
of a telecom fiber. The fiber is linearly birefringent, and 
also characterized by both conservative and dissipative 
cubic nonlinearities. The main feature that makes our 
theory different from most previous studies on fiber-optic 
Raman amplifiers is its vectorial nature. Thus, we care- 
fully consider the propagation dynamics of two polariza- 
tion components of each of the two beams. In total, the 
number of field components is four, and they all interact 
with each other via cubic nonlinearity. The first vectorial 
theory of Raman effect in randomly birefringent optical 
fibers was developed by Lin and Agrawal in Ref. [l9| 
and applied to the regime of interaction characteristic to 
what we call here "standard Raman amplifiers". Here 
we arc interested in a totally different regime, namely 
the regime of Raman polarizer. The difference between 
the two regimes is explained below, in the beginning of 
section II. 

We start from the equation of motion for the sig- 
nal field, written for the two-component field vector 



reads 



U s = (u 



SXJ 



) T , where u sx and 

iisy are the amplitudes 



of the normal polarization modes and e y of the fiber: 
U, — Usx&x ~t~ u S yGy. This equation is derived under the 
(as usual for nonlinear optics) unidirectional and slowly 
varying approximations, see for instance (l9l [20j . and 



id z U s + ip'{u) p )d t U s + AB(u s )U s 



3 y s 



U S )U S 



U p )U t 



u s )ui 



+ (u p ■ u s )u; + (u s ■ u;)u p ] 



+ie s g(U; ■ U S )U P = . 



(1) 



A similar equation (with indices p and s interchanged) 
arises for the pump beam, which is characterized by 
the field vector U p . Here 7 SS and j sp are self- and 
cross-modulation coefficients, whose values depend on 
frequency, and therefore in principle are different for the 
signal and pump beams. They are equal to the frequency- 
dependent Kerr coefficient of the fiber. For simplicity 
we assume 7 SS = 7 pp = 7 ps = 7. /3'(<jj p>a ) is the in- 
verse group velocity of the pump/signal beam. e s = 1, 
e p = — cj s /w p . AB(uj PjS ) is the birefringence tensor. For 
a linearly birefringent fiber it takes the form AB(u> p ^ s ) = 
A/3(w PiS ) (cos#<73 + sin^di), where A/3(w PjS ) is the value 
of birefringence at frequency and 9 the angle of ori- 
entation of the axis of the birefringence with respect to 
the reference frame defined by polarization modes e x and 
e y . 03 and o\ are the usual Pauli matrices. 

The orientation angle 9 is randomly varying in fibers. 
In principle, the magnitude of the birefringence Ad also 
varies stochastically. However, as noticed in Ref. [21(, the 
two approaches, one in which 9 is the only stochastic vari- 
able, and the second, where both 9 and A/3 are stochastic 
variables, produce nearly identical results. Thus, here we 
shall develop our theory by assuming the single stochas- 
tic variable 9. Our theory can be seen as a generalization 
of the one beam linearly birefringent theory of Wai and 
Menyuk from Ref. [2l| to the case of two beams interact- 
ing via the Kerr and Raman nonlinearity in a fiber. The 
angle 9 is driven by a white noise process d z 9 = ge(z), 
where (ge(z)) — and (gg(z)gg(z')) — 2L~ x 5(z — z'). 
Here L c is the correlation length, that characterizes the 
typical distance at which 9 changes randomly. 

Details of the theory are presented in Refs. plHI- Here 
we quickly drive through the major steps of this theory 
and show how to obtain the final result - a set of four cou- 
pled first-order ordinary differential equations, one equa- 
tion for each polarization component of two beams. It 
is instructive to present these four equations as two vec- 
torial equations for the Stokes vectors of the pump and 
signal beams. Each Stokes vector has three components. 

(p) o(p)\ 



Namely, the pump Stokes vector S^ p ^ = (S[ 
has components s[ p) = * pl * P 2 + * P i* 
?0») - lift . 12 



jO) 



-pi 



P2 - 
| S (P)| 



p2 



*:i*«a - = l*„il 2 - If„ 2 l 2 , and power 



Similar expressions define the signal Stokes 



vector S( s > . Here the field vector ^ is related to the orig- 
inal Jones field vector U by the relation ^ p . 



J p,s^ P ,s : 



where the 2x2 matrices T v s with elements 
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obey the stochastic differential equations 

±id z T p +AB{u> p )T p = 0, 
id z T s + AB{co s )T s =0, 



(3) 



(4) 
(5) 



where plus (minus) sign stands for the co- (counter-) 
propagating regime of propagation of the two beams, and 



(6) 
(7) 



Here, 9 Z is the derivative of 9 with respect to z. It is 
different from zero owing to the random changes of ori- 
entation of the birefringence axes. Now the polarization 
components of each beam appear to be defined with re- 
spect to the local birefringence axes, while these axes 
rotate stochastically along the fiber length driven by the 
noise source gs(z). 

These transformations eliminate the birefringence 
terms from the equations of motion of ^ p and *f/ s and 
bring about a vast number of cubic terms composed of 
different combinations of ^ p i, ^ p2 , ^si, ^s2 and their 
complex conjugates. Factors in front of these terms are 
products of two coefficients of the form u m u n , or u* m u n , 
or mJ,!i* , where m, n = 1, . . . , 14. Products with m = n 
we shall call self-products, while with m n cross- 
products. Here, m — | ai | 2 — |a2 1 2 , u 2 = ~{a\a 2 + aTa 2 ), 

*2 
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u 3 = i{a\a 2 — o,\a 2 ), u 4 ~ 2a\a 2 , 115 = af 
u 6 = -i(a\ + a 2 2 ), u 7 = a\bi — a 2 b 2 , u s = — (bia 2 + b 



a%bl, u 13 



i(a*bi + a 2 b 2 ), t»u = 
«(oi&i + a 2 b 2 ), 



u 9 = i(bia 2 - am), u w 
ai&2 + bia 2 , u 12 = aibi 
M14 = i(a\b 2 — a 2 bi)- 

In the thus obtained equations of motion for ty p and *5f s 
we perform the ensemble average (over different realiza- 
tions of the random process which describes linear bire- 
fringence). Thus, we write {u m u n ) instead of u m u n . This 
change holds true only in the limit when the stochastic 
variations are faster than the nonlinear beam evolution. 
This is exactly the place in the derivation where our sin- 
gle approximation comes into play. At this point we also 
need to apply the ergodic theorem 



(/) = lim - f dz' f(z'). 

z-hx> Z Jo 



(8) 



Our goal is to calculate ensemble averages of all necessary 
self- and cross-products: in this way we may complete the 
derivation of the differential equations for *S? p and *S? S . 

The equations of motion for u n with n = 1, . . . , 14 
can be easily formulated basing ourselves on equations 
dl]) and ([S]). As these equations are linear, in order to 
find an ensemble average of any function of these coeffi- 
cients (in our case pair products) we need to construct a 
generator. We refer to the Appendix in Ref. 2l| for de- 



With this generator we are able to formulate the equa- 
tions of motion for the ensemble averages of the products 
of the coefficients. Thus the solutions to the equations of 
motion 

3 Z G X = -2L-\G X -G 2 ), (9) 
d z G 2 = 2L-\Gx - G 2 ) t 4A/?K)G 4 , (10) 
3 Z G 3 = ±4A/?(^)G 4 , (11) 
8 Z G A = -L-'G A ± 2A/3(oj p )(G 2 - G 3 ) (12) 

yield the result for the self-products {(uf), (u 2 ), (u 3 )}, 
{(Re 2 M), (Re 2 K)), (Re 2 (u 6 ))}, 

and {(Im (u 4 )), (Im (1/5)), (Im (u 6 ))}, if we associate 
them with {Gi, G 2 , G3} with initial conditions given as 
(1, 0, 0), (0, 1, 0), and (0, 0, 1), respectively. 

The remaining self-products 

{(Re 2 M), (Re 2 M), (Re 2 ( Ug ), (Re 2 ( Ul0 ))}, 
{(Im 2 K)), (Im'M), (Im 2 ( U9 ), (Im 2 ( Ul0 )}}, 
{(Re 2 ( Ull )), (Re 2 (« 12 )), (Re 2 (u 13 ), (Re 2 ( Ml4 ))}, 
and {(Im 2 (Mn)), (Im 2 (ui 2 )), (Im (u 13 ), (Im (u u ))}, 
can be found from the equations 

d z G x = -2i- 1 (G 1 -G 2 )+2A±G 5 , (13) 

d z G 2 =2L- 1 (G 1 -G 2 )-2A±G 6 , (14) 

d z G 3 = 2A±G 6 , (15) 

d z G A = -2A±G 5 , (16) 

d z G 5 = A±(G 4 -G 1 )-i- 1 G 5 , (17) 

9 2 G 6 = A±(G 2 -G 3 )-L- 1 G 6 , (18) 

when we associate them with {Gi, G 2 , G3, G 4 }, 
with initial conditions as (1, 0, 0, 0), (0, 0, 0, 1), 
(0, 1, 0, 0), and (0, 0, 1, 0), respectively. Here A± = 
[±A/3(uj p ) - A/3(u) s )], 

In order to find the cross-products we constructed 
appropriate generators and found that all the cross- 
products that are of interest to us turn out to be equal 
to zero. Similarly, terms of the form Ke(u n )Im(u n ) also 
vanish. Thus, many SPM, XPM, and Raman terms in 
the final equations of motion disappear. The thus found 
equations of motion for the fields are conveniently formu- 
lated in Stokes space. They read as 

(±a z + /3'K)9 4 )s (p) = 

7 (sto x J^)(z)S (p ) + x J x (z)S( s A 

+e p (g/2) (s^JroS^ + Jr(z)S^) , (19) 

(9 z +/3'K)9 t )S (s) = 
7 (s^ x 4 s \z)S^ + x J x (z)S^ 



tails of this procedure, and only give here the final result. 



+(g/2) (s^J RQ S^ + sPMz)SW) . (20) 

Matrices in equations (|19p and ([2U)l are all diagonal 
with elements 

J R = diag(jRi, J R2 , J R3 ), J x = diag(J xi , J x2l J x3 ), 
J s = diag(J s i, J s2 , J s3 ). These elements are different 
for the counter-propagating and the co-propagating in- 
teraction geometries. 
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In order to complete our theory, we need to express all 
elements in these matrices in terms of ensemble averages 
of self-products: 



Jri — 

JR3 — 


(Re(u 2 4 - u? )) , 
— CRe(v?< 4- w?„^ 

\ 1 >-'-V u 14 • a W)l J 
-(|«14| 2 -Ko| 2 ), 




(21) 
(22) 
(23) 


Jxl = 


2 

-(Re(u 10 + u 13 -'U9 - 




(24) 


Jx2 = 


2 

-(Re(u 10 + u 14 -u 9 - 


«? a )>, 


(25) 


Jx3 = 


|<|u 9 | 2 + |wi 4 | 2 - |wi 3 


2 -Kol 2 ), 


(26) 


Jsl — 






(27) 


Js2 — 






(28) 


Js3 = 


I m) i] , 




(29) 



and also J^o = (l u io| 2 + |i«i4| 2 )- Note that our model 
reduces to the one-beam theory of Wai and Menyuk when 
the coefficients uj through U14 are set to zero. 

The Stokes representation is particularly appealing in 
the context of the problem that we are considering. As 
we are interested in the polarization properties of the 
outcoming signal beam, the Stokes vector quite clearly 
presents the polarization vector on the Poincare sphere. 
The evolution of the Stokes vector draws a trajectory of 
its tip on the sphere. Another quantity of interest is the 
degree of polarization (DOP). In those cases where we are 
dealing with an ensemble of beams, the DOP character- 
izes the length of the average Stokes vector. Here again 
the Stokes representation appears to be rather useful. 

Thus, the equation of motion for the Stokes vector of 
a CW signal beam is 

d z S (s) = 7 S (s) x J s (z)S (s) + 7 S (s) x J,(z)S (p) 

+ e p (g/2) [s^S^ + S« J R (z)S^] . (30) 

(With Jrq = 1, which is the case for all situations con- 
sidered below.) Here J s is the self-polarization modu- 
lation (SPolM) tensor, J x - cross-polarization (XPolM) 
modulation tensor, Jr - Raman tensor. All they are di- 
agonal. Elements of these tensors are dependent on the 
magnitude of the birefringence both at signal and pump 
carrier frequencies, that is on the beat lengths Lb{oj s ) 
and Lb(oj p ), and also on the correlation length L c . All 
these three lengths do not exceed 100 m in conventional 
telecom fibers. The physical meaning of each tensor fol- 
lows from its definition. Thus, the SPolM tensor defines 
how two polarization components belonging to the same 
beam interact in the Kerr medium, and thereby rotate 
the Stokes vector. The XPolM tensor has similar mean- 
ing, but now the rotation is due to the interaction of 
polarization components belonging to different beams. 
Finally, the Raman tensor defines polarization-sensitive 
amplification of amplifier. This tensor is of particular 



importance to us. For instance, when all elements of this 
tensor vanish, the Raman amplifier becomes insensitive 
to the SOP of the pump beam, so that we are dealing 
essentially with a scalar model. Conversely, when the 
diagonal elements of the Raman tensor have appreciable 
values, then the theory must be necessary vectorial. 

Certainly, the evolution of the Stokes vector sensitively 
depends on how the elements of these tensors evolve with 
distance. In order to find their dynamics it is necessary 
to solve the set of linear ordinary differential equations 
which is given above, see also Refs. Instead of writ- 

ing them down here, we present their approximate ana- 
lytical solutions. Fig. QJa,b,c) shows how well these ana- 
lytical solutions reproduce the exact situation. Fig. QTa) 
shows that the elements of the SPolM tensor drop very 
fast and already vanish within the first 10 m of the fiber. 
Given, that the length of the Raman amplifier exceeds 
1 -f- 2 km, we can safely set 



J s = diag(0, 0, 0) . 



(31) 



The elements of the other two tensors also deceases with 
distance, however much slower, namely as 
g 

J x = --diag(l, 1, l)exp(- z/Ld) , (32) 

J R = diag(l, 1, 1) exp(-z/L d ) . (33) 

As demonstrated in Fig. [TJb,c) the decay distance is in- 
deed determined by the characteristic length which 
is called the PMD diffusion length: L^ 1 = ±(D p Aoj) 2 , 
where D p = 2\f2iT\fL~ c / [Lboj p ) is the PMD coefficient, 
[2lj |. and throughout the paper Alu = uj p — uj s is taken 
to be equal to the Raman shift Awjj in the germanium- 
doped silica fibers, that is 13.2 THz. The theory that we 
are developing here is strictly valid only in two limits - 
the limit which we call here Manakov limit [L nl , Lr *C 
Ld) and diffusion limit (Lnl, Lr Ld), where Lnl 
is nonlinear length, and Lr characteristic amplification 
length. 



III. RAMAN AMPLIFIERS VERSUS RAMAN 
POLARIZERS 

Raman amplifiers, which we call here standard Ra- 
man amplifiers, operate in the diffusion limit, as they are 
based on fibers with large PMD coefficients. Thus, for 
D p = 0.2 ps/Vk^ and Aw = Aujr = 13.2 THz, the PMD 
diffusion length Ld is less than 10 m. Taking into account 
that standard Raman amplifiers arc 10 or more kilome- 
ters long, the contribution of the polarization-dependent 
gain (second term in brackets in Eq. (|30)) ) to the total 
gain (both terms in brackets in Eq. (|30|) taken together) 
is totally negligible. The model equation for the signal 
beam is then 



d z S^ = (g/2)S^S 



(34) 



Thus, each component of the Stokes vector is amplified 
independently and equally with the other components. 
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of fast stochastic polarization dynamics of each Stokes 
vector. Quite to the contrary, the Raman contribution 
appears exactly as in the case of isotropic fibers (i.e. in 
absence of the birefringence, and its stochasticity), be- 
cause the mutual polarization scrambling of the relative 
orientations of the pump and Stokes vectors is very in- 
efficient when the PMD diffusion length Ld is long. In 
other words, Raman amplification is insensitive to the 
absolute orientation of the individual SOPs of the signal 
and pump beams in the laboratory frame. It is however 
sensitive to their mutual orientation. In the case of stan- 
dard Raman amplifiers, the signal Stokes vector rotates 
rapidly around the pump Stokes vector, and therefore on 
average "feels" no polarization dependence. In the case 
of Raman polarizers, still the two vectors stochastically 
rotate in the laboratory frame, but they do it now in uni- 
son, so that their mutual orientation is almost "frozen" . 



IV. AN IDEAL RAMAN POLARIZER 
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FIG. 1. Three diagonal elements of the (a) SPolM, (b) 
XPolM, and (c) Raman tensors. In figures (b,c) all three 
curves visually coincide; the blue curve is the analytical re- 
sult showing the exponential decay: oc exp(— z/Lj). Param- 
eters are: L c = 1 m, Ls(oj s ) = 10 m, uj p — uj s — 13.2 THz, 
\ s — 1.55 /im, and A p = 1.45 /im. The PMD diffusion length 
is Ld — 870 m. Note that a brief transient in Fig. \V[b) is not 
resolved on the chosen scale. 



For such Raman amplifier there is no preferentially am- 
plified polarization mode. The model is essentially a 
scalar one. 

A different situation arises in the Manakov limit. For 
PMD coefficients less than 0.02 ps/Vkm, the PMD dif- 
fusion length becomes greater than 1 km. In this case we 
can write the model equation for the signal Stokes vector 
in the form 



<9 z S (s) = - 7 S^ x S< p > 
+ (.9/2 



S (p) S (s) 



S W S (P) 



(35) 



with 7 = I7. In this limit (Ld — > 00) we deal with 
an ideal Raman polarizer. The equation above includes 
two contributions. The XPolM contribution is a cross- 
phase modulation (XPM) part of the Manakov equation, 
in which the factor of § appears as the result of averaging 



As characteristic to isotropic fibers, the signal expe- 
riences maximal gain when its Stokes vector is aligned 
along the pump Stokes vector. To show this we can 
choose (without loss of generality) the pump Stokes 
vector be aligned along its first component: S^ p ^ = 
Sq P \i, 0, 0). Then, we may write for the signal first 
Stokes component: 



d z s[ e) = (g/2)S^ 



S, 



(36) 



If initially the signal Stokes vector is also aligned with its 
first component, then the signal amplification coefficient 
is g. This value should be contrasted to the value of g/2, 
which is characteristic to standard Raman amplifiers, see 
Eq. dMD- 

The other two components of the signal Stokes vector 
are amplified less efficiently than the first component. 
Indeed, their equations of motion are: 



d z S. 



(p) oW 



(g/2)S^S. 



2 



(g/2)S^si s) 



(37) 
(38) 



Here, the gain is only g/2. 

The observations derived from Eqs. (|3"6"1) - (|3"8"|) explain 
the ability of a Raman polarizer to re-polarize light. They 
demonstrate that only the Stokes component of the sig- 
nal aligned with the pump Stokes vector is dominantly 
amplified. In a high-gain Raman amplifier, the differ- 
ence in gain for polarization components may become 
so large that the polarization of the outcoming beam is 
almost perfectly aligned with the pump SOP. This ef- 
fect of alignment is called polarization attraction, or po- 
larization trapping. Shortly, we shall quantify effect of 
the polarization attraction in terms of the DOP, the so- 
called alignment parameter, and some other parameters, 
while now we comment on the output SOP of the outcom- 
ing signal beam measured with respect of the laboratory 
frame. 
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As we have seen, the Raman tensor decays as the dis- 
tance grows larger, see Eq. (1331) . Therefore, it is prefer- 
able to decrease the total fiber length at the price of 
increasing the pump power. Indeed, the first proof-of- 
principle experiment reported in Ref. [16j by Martinelli 
et. at, was carried out with a dispersion-shifted fiber 
of only 2, 1 km, and an average pump power as high as 
2.2 W. 

Most theories of Raman polarizers reported so far, see 
Refs. HiHli, are based on numerical simu- 
lations of the stochastic equations which properly take 
into account the randomness of the fiber birefringence 
(the total fiber span is divided into short segments, with 
each segment extended over one correlation length; the 
orientation of the birefringence axes is fixed within each 
segment, while it varies randomly when going from one 
segment to the next one). Such an approach is indeed 
necessary when the PMD diffusion length is comparable 
with the amplification length, a case which is in the mid- 
dle between the Manakov limit and the diffusion limit. 
In this case, the Raman polarizer has non-optimal per- 
formances, yielding a DOP which is significantly below 
unity. So, this case is not advantageous in practice. In 
order to realize a "good" Raman polarizer, one should 
choose to work in the Manakov limit. As we have indi- 
cated above, working in this limit allows us to get an- 
alytical and physically transparent results. In the next 
section we shall continue to work with ideal Raman po- 
larizers and provide an even deeper analytical insight. 



V. EVALUATION OF THE PERFORMANCE OF 
A RAMAN POLARIZER 

Equations (|3"o]) - (13"51) can be solved analytically. We 
shall limit ourselves to the undepleted pump approxi- 
mation, so that the pump power P = Sq(z) =const. 
Our model does not include linear losses in the fiber, be- 
cause we have chosen to work with relatively short fiber 
spans, for which losses are relatively small. If necessary, 
the losses can be included, though analytics will become 
less transparent. Solutions to Eqs. (|3"o]) - (15"5|) are: 



si s \z) = ~[si s \o)-s[ s \o) 



1 

+ 2 



S ( o s) (Q) + S i 1 s) (0) 



S<j s) (0) + S< s) (0) 



si s \z)=[si s \0)co S (jPz) 
-S ( 3 S) (Q) sin(jPz) 



\gPz 



+S ( 3 \o) cos^Pz)]ei 9Pz . 



(39) 



(40) 



(41) 



(42) 



We are interested in the statistical properties of a Ra- 
man polarizer. All quantities of interest can be derived 
from the above-written solutions. First of all, we shall 
calculate the mean quantities. The immediate questions 
are - what is the SOP of the outcoming signal beam and 
how well the beam is polarized? In order to find an an- 
swer to the first question we need to simply get an aver- 
age of Eqs. (13T)1) - (|4"2")) given the statistics of the incoming 
light. We assume that the signal is initially unpolarized, 
so that (S^ s) (0)) = (S { 2 S) (0)) = (S { 3 \o)) = 0. Then, at 
z = L, where L is the total length of the fiber, we get 



(S^(L)) = ^ s \0) [1 + expfoPZ)], 

{4'HL)) = \sP(p) [-1 + expbPL)] , 

(S [ 2 s \L))=0, 
(S { 3 s) (L))=0. 



(43) 

(44) 

(45) 
(46) 



So, the signal SOP at the output is aligned with the 
pump SOP. The degree of alignment is characterized by 
the DOP, which is calculated as 



J (S[ s) (z)> 2 + (S ( 2 S) (z)) 2 + (S { 3 S) (z)) 2 

DOP(z) = V , 2 . 3 . (47) 

(S ( s) (z)) 

As usual, a DOP equal to unity means that light if 
perfectly polarized, a DOP equal to zero indicates that 
the light beam is unpolarized, while intermediate val- 
ues stand for a partially polarized beam. We say that 
the Raman polarizer perfoms its function properly when 
DOP becomes close to unity. Introducing gain G as 
G = (S^ s \L))/S { o s \0) we get G = ± [l + exp(gPL)] and 
for the DOP: 



DOP = 1 — G~ 



(48) 



The higher the gain, the larger the DOP. Already 20 dB 
gain is enough to get a DOP as high as 0.99. 

A short comment is in order on how one should 
interpret the averaging procedure, expressed by (...). 
There are two possible situations. On the one hand, 
we can vary the SOP of the signal beam in time, then 
(...) = (. . .) T = T- 1 ... dt, where T is the period 
of time, sufficiently long to get correct statistical averag- 
ing. (s[ s ^)t = (i?2 ) T ~ (S^^t — means that we are 
dealing with unpolarized light. On the other hand, we 
can imagine an experiment with an ensemble of beams. 
Then, (...) = (. . .) e means ensemble average over all 
these beams. If the SOPs of all beams from the ensemble 
randomly or uniformly cover the Poincare sphere, then, 
similarly to the time average, we get (S[ s ^) e — (S^'^e — 
(S 3 s ^} e = 0. In this situation we say that we are dealing 



with an ensemble of scrambled beams. In a case where 
the time average gives the same statistical information 
as the ensemble average, we refer to such system as an 
ergodic one. The Raman polarizers considered here are 
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obviously ergodic systems, simply because time does not 
enter the equations of motion explicitely. Therefore, our 
analysis is valid for the scrambled beams approach as well 
as for time averaging. 

Another important quantity which characterizes a Ra- 
man polarizer is the alignment parameter A^f, defined as 
the cosine of the angle between the output signal SOP 
and the output pump SOP: 



.4 



tt 



a( s ) qip) 
D 2 °2 



(49) 



The closer the alignment parameter to unity, the bet- 
ter the alignment of the output signal and pump Stokes 
vectors. Using solutions in Eqs. (|43l) - (|4l))) we get 



(s[ s \l)) 

(S ( 



= 1-G' 



(50) 



for the value of the alignment parameter at the fiber out- 
put. Although this value coincides with the value of 
DOP, see Eq. (|48|). these two quantities have different 
physical meanings. For a statistical ensemble of beams, 
the alignment parameter shows the average direction of 
the signal Stokes vector on the Poincare sphere, while the 
DOP measures the breadth of the spot traced by the tips 
of the signal Stokes vectors around this average direction. 

Yet another quantity of interest is the measure of 
the polarization-dependent gain (PDG). It is exactly the 
PDG which is at the heart of a Raman polarizer. Dif- 
ferent SOPs of the signal beam experience different am- 
plifications. The signal beam with a SOP parallel to the 
pump Stokes vector is amplified most efficiently, while the 
orthogonal polarization experiences no gain. Indeed, as it 
follows from the solution in Eq. (|40p. G max = exp(gPL) 
and G m in = 1- We introduce the PDG parameter A as 
A = G max — G m ; n , and get for the ideal Raman polar- 
izer A = 2(G — 1). The PDG parameter aquires high 
values for a high-gain Raman polarizer. Note that for 
an "ideal Raman amplifier" (an amplifier, which is per- 
fectly described by the scalar theory, or in other words, 
the amplifier, which works deeply in the diffusion limit) 
A = 0. 

The high value of the PDG parameter points out 
that along with the desirable property of strong re- 
polarization of the signal beam, this device is character- 
ized by a high level of unwanted relative intensity noise 
(RIN). By varying the signal SOP at the input we get pro- 
nounced variations of the intensity at the output, even 
if the incoming beam had a steady intensity in time. In 
order to characterize the output power fluctuations, let 
us calculate the variance 



1 



Using solution in Eq. (|39l) we get 

a? = (l-G- 1 ) 2 /3 



(51) 



(52) 



For large G, a s « « 58%.This level of RIN may be 

detrimental for some optical devices, particularly nonlin- 
ear ones. Note that an ideal Raman amplifier is charac- 
terized by a s = 0, i.e., by zero RIN, thanks to the efficient 
polarization scrambling which is provided by PMD. The 
price to be paid is the totally stochastic signal SOP at 
the output fiber end. 

A reasonable question to ask is whether it is possi- 
ble for a Raman polarizer to keep the useful property 
of re-polarization and at the same time to suppress RIN 
down to an acceptable level. The answer is positive. One 
possible way to combat the RIN and at the same time 
keep the property of re-polarization is to use the Ra- 
man polarizer in the depleted-pump regime, Ref. [24j . In 
this saturation regime all input SOPs are amplified to 
approximately the same level of intensity, actually up to 
Sq S \l) P. Strictly speaking, only one signal SOP (the 
one which is perfectly orthogonal to the pump SOP) is 
not amplified at all. However, the numerous imperfec- 
tions of any practical realization of a Raman polarizer, 
including residual PMD, may prevent the observation of 
such a singular behaviour. 

So far we have analyzed the main statistical properties 
of an ideal Raman polarizer operating in the undepleted- 
pump regime. If necessary, any other statistical quantity 
of interest can be obtained from the exact analytical so- 
lutions given in Eqs. (|39p - (|42[) . In a similar manner, one 
can characterize the re-polarization of partially polarized 
beams. The final quantity which we would like to com- 
ment on is the mean gain of an ideal Raman polarizer. 
It is well known, that the gain of an ideal Raman ampli- 
fier is equal to g/2. The reason is that in the course of 
propagation the signal SOP rotates quickly around the 
pump SOP, and on average "feels" the arithmetic mean 
of the maximal gain (g) when it is parallel to the pump 
SOP, and minimal gain (0) when it is orthogonal, yield- 
ing g/2 on average. In terms of available gain, an ideal 
Raman polarizer performs much better. As can be seen 
from Eq. (|43|) . for large values of G, G sa exp(gPL — ln2), 
so that the gain coefficient is almost twice larger. This 
property makes Raman polarizers very efficient Raman 
amplifiers as well. 



VI. COUNTER-PROPAGATING RAMAN 
POLARIZERS 

So far, we have been dealing only with the co- 
propagating geometry. In this geometry, the pump SOP 
stochastically changes along the fiber, and its output 
SOP depends on the particular realization of the bire- 
fringence stochasticity in the chosen fiber span. More- 
over, the stochasticity changes with time, as a result 
of variation of the environmental conditions. Therefore 
the trapping of signal's SOP to pump's SOP does not 
garantee the absence of fluctuations of signal's SOP at 
the output, even though these fluctuatons closely fol- 
low the time-varying pump SOP. In other words, the 
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co-propagating Raman polarizer provides the trapping 
effect in the stochastic frame, but does not garantcc the 
SOP stabilization in the laborotary frame. 

The desirable stabilization in the laboratory frame can 
be achieved by implementing a counter-propagating ge- 
ometry, Refs. @, |23|- Since the signal's SOP is attracted 
toward the instantaneous position of pump's Stokes vec- 
tor, this alignment holds also at the output end of the 
fiber. The output pump SOP is defined solely by the 
source, and as such it is supposed to be well defined and 
deterministic. In this respect the counter-propagating ge- 
ometry is prefcrrable. As regarding the theory, one can 
repeat derivations with the opposite sign of z-derivative 
in the equation governing evolution of the pump beam. 
As shown in Ref. Q, this reversing of the sign brings 
some changes in the components of the XPolM and Ra- 
man tensors. They become 

jcounter = _^ diag{h 1} CX p(-z/L d ) , (53) 

jcounter = _^ y exp („ z/L(i ) . ( 54 ) 

The presence of the factor | in front of the Raman tensor 
immediately leads us to the conclusion that the counter- 
propagating Raman polarizer is significantly less effective 
in re-polarization than its co-propagating analog. In or- 
der to get similar performances we need either to increase 
the pump power or lengthen the fiber, or both. Let us 
evaluate the performance of this device. 

First of all, we start with the solving the equation of 
motion (|35l) in the undepleted-pump regime. We get 

1 r 



l r 
2 



S ( s) {0)- S\ s >(0) 



(«), 



s { s) (0) + s[ s) (0) 



1 r 



1 r 
2 



S ( s) (0)- s[ s \o) 



\ a Pz 



s[ s) (0) 



ef " p * . 



SP(z)= SP(0)cos(jPz) 



-S { 3 s) (0) sm(jPz) 



s£°(0) sin( 7 Pz) 



4 s) (0)cos(7Pz) 



C2 



(55) 



(56) 



(57) 



(58) 



We immediately observe that the difference in amplifi- 
cation coefficients of the first Stokes component and the 
second (and third) Stokes component is given by | g— ^g. 
The contrast is much weaker than for the co-propagating 
case, where we had g — \g. The average gain of the 
counter-propagating Raman polarizer is 



G = 



C 39PL + giffPi 



(59) 



which is significantly smaller than for a Raman polarizer 
operating in the co-propagating configuration, although 



it is still larger than for an ideal Raman amplifier. For the 
same value of the product PL, the DOP for the conter- 
propagating configuration is also smaller: 



DOP = 1-2 



1 



1 - 2e~i gPL [iovgPL > 1) 
1 - V2G- 1 / 2 . 



(60) 



For G = 20 dB in the co-propagating case the DOP was 
as high as 99%, while in the counter-propagating config- 
uration it is only 86%. 

It is instructive to compare our model of ideal Raman 
polarizer with full-scale numerical simulations of the un- 
derlying stochastic equations presented in {23} , where the 
empirical formula: 



DOP = 1 - e ~ GdB/T . 



(61) 



connecting the DOP with the gain was suggested and 
tested numerically. Here GdB = 101og 10 G and Y w 10.2 
for the considered range of PMD coefficients. The graph- 
ical comparison of the results obtained with formula (|60[) 
on one hand, and the results plotted according to the 
empirical formula (|61[) on the other hand, is shown in 
Fig-UJa). The fit is good. On the same plot we have also 
shown the results based on the direct numerical solution 
of Eq. (|3"0")) with XPolM and Raman tensors in the form 
of Eqs. (|5"3"|) - (l5^i|) . Note that we did not use any fitting 
parameter in this cross-comparison. 

The alignment parameter for the counter-propagating 
geometry is different from the co-propagating case. Be- 
cause of the change of the sign in front of the second el- 
ement of the Raman tensor, see Eq. (15^]) . the alignment 
parameter is now 



.4 



(s{ s) s[ p) 



c( s ) q(p) 
°2 °2 



(62) 



For input unpolarized light, the alignment parameter co- 
incides with the DOP, namely, 



A n = l-2[e* 



i 



The PDG parameter A = G 
lated, resulting in 



1 - 2c^ 9Pi (ioigPL > 1) 

1 - V2G~ 1/2 . (63) 

max - G mi „ is easily calcu- 



A 



1 / 2 

- ea 

2 V 



aPL _ e igPL 



1 



1 + 2G-V1 + 8G 



(64) 

Its value is considerably less in the co-propagating con- 
figuration, particularly for moderate values of gain. This 
observation again points to the relatively poorer perfor- 
mances of the counter-propagating Raman polarizer. At 
the same time, the RIN is expected to have a lower level. 
In order to demonstrate this, let us evaluate the variance 
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53T) . and (l54l were calculated in the limit 
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FIG. 2. DOP versus gain G. Graphical comparison of the 
results obtained with formulae (|60[) (black solid); results ob- 
tained with empirical formulae (|6ip (green dotted); and re- 
sults based on the direct numerical solution of Eq. (I30|l with 
XPolM and Raman tensors in the form of Eqs. (I53p ~ (|54p (red 
dashed). Parameters are: a) L c — 1 m, Lb — 45 m, P = 8 W, 
D p = 0.005 ps/x/km, and Ld = 17.5 km, L varies from km 
to 2.5 km, and T = 10.2; b) L c = 10 m, L B = 3500 m, 
P = 8 W, D p = 0.0002 ps/Vkm, and L d = 10914 km, L 
varies from km to 1.5 km, and T = 4.3. 



of the signal intensity, 
yield 



Formula (1ST1) and solution ([55 



J 1 



1 - 2 



gPL 



1 2 



(65) 



Before concluding this section, one remark is in order 
about the applicability domain of these results. SPolM, 
XPolM and Raman tensors given by Eqs. (HP, (J32]), (l33l) . 



L > L b 



T2 

B 



(uj p ) 



8ir 2 L c 



(66) 



This inequality holds for all practical situations. Thus, 
for L c as small as 1 m and Lb as large as 100 m we get 
Lfcire a $ short as 127 m. Fiber-optic Raman amplifiers are 
always longer than 1 km, and therefore inequality (|66|) is 
not violated. However, if for some reason inequality (|66[) 
is violated, for instance for extremely low birefringent 
fibers, the analysis given above must be corrected. Thus, 
in the limit Lb(uj p ) — > oo, the tensors of interest take the 
following form: 



jcounter = di&g (_ 
jcounter = 4^ 



■r 



R 



1, 1, -1), 

~1,0, 1), 
diag(l, 1, 1) . 



(67) 
(68) 

(69) 



Fig- Hlb) shows the dependence of DOP on Raman polar- 
izer gain for this case. Although the performance of the 
Raman polarizer in this limit is very good, we will not 
evaluate it here explicitely because of its little practical 
interest. 



VII. CONCLUSION 

We have studied the effect of trapping of the state 
of polarization of a signal beam by a pump beam in 
the model of a Raman polarizer. We have introduced 
the notion of the ideal Raman polarizer and quantified 
its performance it terms of gain, degree of polarization, 
polarization-dependent gain parameter, alignment pa- 
rameter, and RIN characteristics. We have studied two 
different geometries: co - and counter-propagating con- 
figurations, and identified their pros and contras. Possi- 
ble applications of Raman polarizers include their poten- 
tial use in telecom-related signal processing, where the 
need of transforming an unpolarized light to a polar- 
ized one is necessary in order to provide an interface be- 
tween the telecom link and post-processing polarization- 
sensitive devices (based, for instance, on nonlinear crys- 
tals) . 
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